Quantum and thermal phase transitions in the Bechgaard salts and their analogs 
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In order to investigate quantum and thermal phase transitions in the Bechgaard salts and their 
sulfur analogs, we perform finite-temperature Hartree-Fock calculations in one dimension with par- 
ticular emphasis on the interplay between charge ordering and lattice instability. The coexisting 
charge- and spin-density-wave state as well as its precursor fluctuations in (TMTSF) 2 PFe and the 
lattice tetramerization in (TMTTF) 2 Re04 are well interpreted. 
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The Bechgaard salts (TMTSF) 2 X and their sulfur 
analogs (TMTTF) 2 A, where X = PF 6 , Re0 4 , Br, etc., 
exhibit various types of charge and/or spin ordering 
W and stimulate our interest in quasi-one-dimensional 
correlated electron systems with a three-quarter-filled 
(quarter-filled in terms of holes) 7r band. A vari- 
ety of mechanisms have been proposed in an attempt 
to interpret the variety of instabilities, featuring di- 
mensionality which can be tuned by pressure as 
well as chemically, intrachain dimerization due 
to the anion columns, and competing Coulomb inter- 
actions — on-site, nearest-neighbor, and even 
next-nearest-neighbor repulsions. Recently, introduc- 
ing another interesting viewpoint— the interplay between 
charge ordering and lattice instability, several authors 
[^HH have revealed further exotic density-wave states. 
In particular, Riera and Poilblanc (||| found charge- 
density-wave (CDW)-bond-order-wave (BOW) mixed 
phases of new type, which consist of the superposition 
of lattice dimerization and tetramerization, and pointed 
out their relevance to the spin-Peierls transition re- 
cently observed in TMTTF salts. 

Thus motivated, we explore thermal, as well as quan- 
tum, phase transitions in the Bechgaard-Fabre salts tak- 
ing account of long-range Coulomb and electron-lattice 
interactions. We consider a one-dimensional model at 
quarter filling, 



H. = "Ht + 7~Ci 



Hy + Ti-K ; 
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FIG. 1. Schematic of various spin-charge-lattice instabili- 
ties, where the various circles, segments, and vertical arrows 
qualitatively describe the variation of charge densities, bond 
orders, and spin densities, respectively, while the horizontal 
arrows denote lattice distortion from equilibrium. 

with Cj a creating a hole of spin a on the jth molecular 
site, Uj being the displacement of the jth molecular site 



(2) 



from equilibrium, rij, a — Cj a Cj tCr , and rij — nj,-\ + nj,i- 
Considering the intrinsic dimerization of the system due 
to the anion columns, we introduce alternation to trans- 
fer integrals, nearest-neighbor Coulomb interactions, and 
spring constants as tj = —t\ + a\(uj+i — Uj), Vj = V\, 
and Kj = K x , where A is set to "a" for an odd j and to 
"b" for an even j. We calculate the free energy within 
the Hartree-Fock approximation, which should be dis- 
tinguished from the mean-field treatment of the Hartree 
level [^Pj9|,[l0| and is crucially superior to it with intersite 
Coulomb interactions included. The lattice distortion is 
adiabatically treated and can therefore be described in 
terms of the electron density matrices. We always set 
the smaller ones of the intrinsic transfer integrals and the 
spring constants, and K^, both equal to unity. The 
Hamiltonian (]l|) possesses various density-wave solutions, 
some of which are schematically shown in Fig. ??, where 
(A), (A'), and (B') are accompanied by spin alignment 
and their nonlocal stabilization into a spin density wave 
(SDW) assumes weak interchain interaction. 
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FIG. 2. Ground-state phase diagrams: (a) t a /tb ~ 1-1, 
U/t h = 4.0, V h /V* = 0.8, a a /(f b ^a) 1/2 = 0.1, 
a h /(t b K b )^ 2 = 0.11, A- a //C b = 1.1; (b) i a /t b = 1.35, 
*7/*b = 4.0, Vb/K = 0.7, a a /(tbif a ) 1/2 = 0.5, 
a h /(t h K b ) 1/2 = 0.6, tf a /tf b = 1.2. 
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FIG. 3. Thermal phase diagrams on the lines of 
V'/Ve. = 0.75 in Fig. 2(a) and of V'/V a = 0.6 in Fig. 
2(b), where the dotted arrows imply the thermal behaviors 
observed in (TMTSF) 2 PF 6 (a) and in (TMTTF) 2 Re0 4 (b), 
respectively. 



is of second order, whereas transitions between (A) and 
(B) are of first order. Therefore, the 2k F CDW and SDW 
fluctuations may in principle coexist in the intermediate 
temperature region. It is true that we narrowly find such 
a multistep thermal behavior within the present model, 
but we could ensure the possibility by including the cru- 
cial effect of the anion potential f§. In (TMTTF) 2 Re0 4 , 
on the other hand, the lattice-tetramerized CDW state 
(B) is stabilized below about 160 K. There is neither 
precursor SDW fluctuation at intermediate temperatures 
nor static spin ordering at low temperatures. Figure ||(b) 
can qualitatively interpret all these observations. Apply- 
ing pressure, or equivalently reducing bond alternation 
by chemical tuning of materials, Chow et al. plq] observed 
crossover from spin-Peierls disordered states to 2/cp SDW 
states in the Bechgaard-Fabre salts. Considering that in 
the 2kp CDW state (B) electrons become more and more 
bound in singlet pairs on next-nearest-neighbor bonds 
with decreasing temperature and therefore the state (B) 
is a fair realization of the observed spin-Peierls phase, 
Figs. ||(a) and ||(b) are widely consistent with experi- 
ments, giving a hint in assigning parameters to the ma- 
terials. 
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We show in Fig. ^| ground-statephase diagrams cal- 
culated for typical parameters §|(J of (TMTSF) 2 PF 6 
(a) and for those @,|| of (TMTTF) 2 Re0 4 (b). The in- 
trinsic dimcrization (bond alternation) is more remark- 
able in TMTTF salts than in TMTSF salts in general. 
Figure |^(a) may be compared with a pioneering calcu- 
lation by Kobayashi et al. j|3], which was carried out 
within the Hartree approximation under frozen lattice. 
The two calculations are qualitatively alike, but the long- 
range Coulomb effect of the Fock field on magnetic phases 
strongly reduce the coexisting CDW and SDW phases 
(A') and (B'). 

Taking account of experimental observations [[[]|l|] , we 
focus on particular regions in Fig. || and show their ther- 
mal behaviors in Fig. § In (TMTSF) 2 PF 6 , a precursor 
CDW fluctuation of type (B) comes at about 150 K, it is 
reduced and a SDW fluctuation of type (A) grows instead 
at about 50 K, and then the mixed state (B') is stabilized 
below 12 K. Figure |^(a) well solves these observations. 
Any transition to (B'), whether from (A) or from (B), 
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